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Abstract 

The  nonparamet r ic  constant  false  alarm  rate  (CFAR)  property 
of  digital  polarity  sampled  phase  reversal  coded  pulse  compressors 
is  described.  The  detection  performance  in  Gaussian  and  non- 
Gaussian  noise  is  determined.  It  is  shown  that  the  loss  in  signal- 
to-noise  ratio  of  the  processor,  relative  to  the  incoherent 
matched  filter,  decreases  as  the  code  length  increases,  the 
asymptotic  loss  in  Gaussian  noise  is  1.96dB,  and  the  loss  in 
Weibull  noise  decreases  with  the  shape  parameter  of  the  Weibull 
distribution  and  can  even  become  a  gain. 

I.  Introduction 

While  doing  pulse  compressing,  digital  polarity  sampled  phase 
reversal  coded  pulse  compressors ^  also  completes  constant  false 
alarm  rate  (CFAR)  processing.  This  type  of  processors  are  easy  to 
Implement  and  thus  are  widely  used  in  radars.  Ref  [2]  has  analyzed 
the  loss  mechanism  of  this  type  of  processors  and  given  the 
quantitative  study  of  the  loss  in  typical  applications  with 
Gaussian  noise  and  finite  code  length.  This  paper  studies  the  non- 
parametric  CFAR  characteristics  and  the  detection  performance  in 
Gaussian  and  non-Gaussian  noise  of  this  type  of  processors.  We 
derive  the  asymptotic  relative  efficiency  (ARE)  in  Gaussian  noise, 
and  determine  the  los3  of  signal  to  noise  ratio  by  using  analytic 
and  numerical  calculations.  Part  of  the  material  has  been 
discussed  in  Ref  [3] . 

II.  Non-parametric  CFAR  characteristics 

The  schematic  diagram  of  the  processors  is  shown  in  Fig.  1. 
When  the  signal  is  present,  the  complex  variable  representation  of 
the  mid-frequency  received  wave  form  is 


2 


(1) 


Zi5*  Vi  e^1  *  A  +•  Ni  e J*t,  iml,  *  •  ■ ,  M, 

where  Vi,  A  and  Ni  represent  received  wave  form,  signal  and  noise 
amplitude  corresponding  to  the  i-th  code  element;  0i,  0i,  and  yi 

indicate  their  corresponding  phases.  Among  them  the  phase  of  the 
signal  8i=0o+cti,  and  0q  is  the  initial  phase.  Depending  on  the 
coding,  <Xi  takes  values  of  either  0  or  it.  M  is  the  code  length. 
The  synchronous  phase  of  Zi  and  the  orthogonal  components  of 
polarity  samples  sgncosai  and  sgnsin<{)i  are  separately  transferred 

to  the  duel  signal  registers  and  are  compared  digit  by  digit  with 
the  codes  stored  in  the  coding  registers.  If  they  match,  +1  is 
chosen.  Otherwise  -1  is  chosen.  Sum  up  the  results  of  the 

comparison  to  get  the  total  of  each  circuit  branches 

M 

I  =  X  sgncos<J>i  sgncosai,  (2) 

i=l 
M 

Q  =  X  sgncosai  sgncosai,  (3) 

i=l 

The  detection  statistics 


t 


I2  +  Q2 
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(4) 


is  obtained. 


Fig.  1  Processor  schematic  diagram.  1.  middle  frequency 
signal,  2.  amplitude  limiter,  3.  self  excitation,  4.  filtering,  5. 
polarity  sampling,  6.  signal  register,  7.  comparison  adder,  8. 
coding  register. 


When  the  noise  phase  \j/i  independently  co-distr ibutes  and  the 
distribution  of  i  does  not  depend  on  the  noise  amplitude 
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distribution,  the  false  alarm  probability  of  this  type  of 
processors  is  independent  of  the  noise  amplitude  distribution, 
i.e.  it  has  CFAR  characteristics.  In  radar  applications,  circular 
symmetric  noise  distribution  is  often  encountered.  Its  amplitude 
and  phase  are  mutually  independent  and  the  phase  distribute 
uniformly  on  {0,  2%)  .  The  sampling  circuit  generally  uses  code 

sampling  and  the  samples  are  independent.  Under  these 
circumstances,  the  processors  have  the  non-paramet ric  CFAR 
capability.  This  type  of  processors  can  also  be  regarded  as  a 
simplified  version  of  Dicke-fix  receivers ^ ,  or  a  symbol  detector 
of  narrow  band  signals ^  . 

III.  Asymptotic  performances  in  Gaussian  noise 

Due  to  the  influence  of  the  noise,  4>i  in  general  is  different 
from  9  j_ .  In  Gaussian  noise  and  non-undulating  signal  situation,  the 

Pete  error  probability  of  the  circuit  branch  I,  i.e.  the 
probability  that  cos<|>i  and  cosQi  has  opposite  signs  ist6^ 

qi-  1  -  <t>(V2S  1  COS0O  I  )  ,  (5) 


where 


<D(x)  = 


1  t 

=-  J  exp  ( -— )  du 

2K  -oo 


is  the  Gaussian  distribution  function,  S  the  input  signal-to-noise 
ratio  for  the  processors  shown  in  Fig.  1.  Similarly  the  Pete  error 
probability  of  the  circuit  branch  Q  is 

q2-  1  -  0(V2S  Isin0ol)  ,  (6) 

k  n 

Since  (Xi  is  either  0  or  n,  when  -cj-  £  9q  and  cos9i  and  cosa* 
take  the  same  sign,  the  probability  of  sgncos<|>isgncosai  equaling  to 
+1  is  1-qi,  and  the  -1  probability  is  qi .  Hence  I  obeys  binomial 
distribution,  and  its  probability  is 

/M\ 

Pi  (n  =  -M  +  2N)  =  V  )  <  l~qi)  Nqj.M-N,  N  =  0,1,  •••,  M.  (7) 

N 

K  n  3n 

When  ~  ^  0o  <~2~  and  cos0i  and  cosai  take  opposite  signs,  the 

probability  of  I  is  the  same  as  formula  (7)  except  n  -  M  -  2N. 
Only  I2  will  be  used  in  the  following,  hence  formula  (7)  can  always 
be  used  for  any  value  of  9q.  Similarly,  the  probability  of  Q  is 
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PQ(n  =  -M  +  2N)  -  (  )  (l-q2)Nq2M_Nr  N  »  0,1,  •••,  M.  (3) 

N 

Since  the  two  orthogonal  components  of  the  narrow  band 

Gaussian  noise  are  independent,  variables  I  and  Q  are  independent. 

For  very  large  M  and  very  small  S,  t  converges  to  Rice 

distribution.  Its  probability  density  function  is 

t  t2  +  4M2S/7l  / - 

P  (t )  =M  exP<- - 2M - )  Io(2t-yS/JI  ),  t£0.  (9) 


)  I0  (2tyS/K  )  ,  t£0. 


When  only  noise  is  present,  S=0,  formula  (9)  becomes  Rayleigh 

distribution.  Its  probability  density  function  is 

t  t2 

P0<t)  =  -  exp  (-  t£0.  (10) 

The  relationship  between  the  detection  probability  Pd  and  the  false 
alarm  probability  Pf  is 
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Pd=  Q(^21n— , 


4M  S 


Here  Q(-)  denotes  Marcum  Q  function. 

For  the  Gaussian  noise,  the  best  parameter  detector  for 
detecting  phase  coded  signals  with  unknown  initial  phase  is  the 
non-interfering  matching  filter.  When  the  input  signal  to  noise 
ratio  is  Sj,  and  the  code  length  is  Ml,  the  detection  probability  Pd 
and  the  false  alarm  probability  Pf  of  the  non-interfering  matching 
filter  satisfy  the  following  relation^7): 


Pd=  Q  (  v  2  In* 


(2MLSL)  . 


In  order  to  make  the  Pd  and  Pf  of  the  two  kinds  of  detectors 
equal,  the  following  has  to  be  satisfied 

=  2  MLSL.  (13) 

% 

Hence,  the  asymptotic  relative  efficiency  of  the  processor  shown 
in  Fig,  1  compared  with  the  non-interfering  matching  filter  is 


ARE  „  lim  =  £ 

S=SL->0  M(Pd,Pf,SL)  jt 
Similarly,  the  asymptotic  loss  is 

t  _  lim  —  _  H 

”  SL(Pd/PfiSL)  2' 
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L«(dB) 


1.96. 


(15) 


For  undulating  signals,  (11)  and  (12)  should  be  averaged 
according  to  the  undulating  model.  Then  formula  (13)  will  still  be 
valid  to  the  averaged  signal  to  noise  ratio.  Therefore  ARE  and  L«> 
do  not  change . 

IV.  Finite  code  length  performance 

The  probability  distribution  of  I2  and  Q2  are  respectively 
Pl2<n  =  m)  =  Pi(n-'fm)  +  Pi(n^-'/M),  m>0, 

PX2  (n  =  0)  =  Pi  (n=0)  (16) 

and 

Pq2  (n  =  m)  =  PQ(n=Vm)  +  Pq  ,  m>0, 

PQ2  (n  -  0)  =  Pq (n=0)  .  (17) 

Under  Gaussian  noises,  I2  and  Q2  are  independent  variables. 
Therefore,  the  probability  of  t2  is 

Pt2(n)  -  PX2 (n) *PQ2 (n) ,  (18) 

where  *  denotes  convolution.  From  the  noise  only  Pt2(n)  we  can 
determine  the  threshold  corresponding  to  the  specified  false  alarm 
probability  Pf.  Afterwards,  based  on  this  threshold  and  the  value 
of  Pt2  (n)  for  both  the  signal  and  noise,  we  can  obtain  the 
detection  probability  Pd  for  non-undulating  signal.  Note  that  Pd  is 
related  to  0q  .  The  detection  performance  of  non-interfering 
matching  filter  can  be  obtained  from  Ref  [7]  .  The  calculation 
result  is  summarized  in  Fig.  2.  Shown  in  the  figure  is  the 
relationship  between  the  signal  to  noise  ratio  loss  and  M  for  the 
processor  of  Fig.  1  relative  to  that  of  non-interfering  matching 
filter  under  non-undulating  signal  and  when  0o=45°,  Pd~0-5,  Pf=10-4 

and  1CT6. 

For  non-Gaussian  noise,  the  calculation  assumes  Weibull  noise 
model.  Since  the  processor  in  Fig.  1  has  the  non-parametric  CFAR 
capability,  the  threshold  corresponding  to  a  given  Pf  value  in 
Weibull  noises  is  the  same  as  that  determined  under  Gaussian 
noises  in  the  above.  When  calculating  Pd,  Monte  Carlo  simulation  is 
used  due  to  lack  of  analytic  solutions.  The  detector  performance 
of  the  non-interfering  matching  filter  under  Weibull  noises  is 
also  obtained  by  Monte  Carlo  simulation.  The  calculation  result  is 
summarized  in  Fig.  3.  Shown  in  the  figure  is  the  relationship 


6 


between  the  signal  to  noise  ratio  loss  and  the  Weibull  shape 
factor  a  for  the  processor  of  Fig.  1  relative  to  that  of  non¬ 
interfering  matching  filter  under  non-undulating  signal  and  when 
0o-45°,  Pd-0.5,  Pf-10-4  and  10“6.  The  tail  of  the  Weibull  noise 

envelop  distribution  becomes  longer  as  a  decreases.  The  a- 2 

Weibull  distribution  is  the  narrow  band  Gaussian  noise  situation. 

We  can  see  that  the  loss  of  signal  to  noise  ratio  decreases  as  a 

decreases.  It  even  becomes  a  gain. 
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Fig.  2  Loss  in  Gaussian  noise.  Fig.  3  Loss  in  Weibull  noise. 

1.  loss,  2.  non-undulating  signal. 

V.  Conclusion 

Polarity  sampled  phase  reversal  coded  pulse  compressors  have 
the  nonparamet ric  CFAR  capability  and  are  easily  realizable. 
Compared  to  non-interfering  matching  filters,  the  relative 
detection  performance  of  this  kind  of  processors  improves  as  the 
code  length  and  the  tail  of  the  noise  envelop  distribution 
increase.  In  Gaussian  noises,  their  asymptotic  loss  is  1.96dB. 
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